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Abstract. Experimentally observed grey solitons are analytically extracted from a 
physically viable Gross-Pitaevskii equation. Associated Lieb and Bogoliubov modes 
arc calculated for these class of solitons. It is observed that, these nonlinear excitations 
have strong coupling with the trap at low momenta and hence can be effectively isolated 
from the Bogoliubov sound modes, which responds weakly to harmonic confinement. 
This strong coupling with the trap also makes the grey soliton amenable for control and 
manipulation through both trap modulation and temporal variation of the two-body 
interaction. 
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1. Introduction 

Recent observation of oscillations between stable localized grey solitons in a cigar shaped 
Bose-Einstein condensate (BEC) and vortex rings [1] has led to significant interest 
regarding interaction between the collective matter wave excitations. These coherent 
structures are relevant for atom interferometry, making their control and manipulation 
an area of active research [1, 2, 3, 4]. The nonlinearity of Gross-Pitaevskii (GP) equation, 
describing BEC at the mean-field level, leads to these excitations, which have been 
explored extensively in recent times [5, 6, 7, 8]. The GP equation reduces to the well 
known nonlinear Schrodinger equation (NLSE) for a cigar-shaped BEC, which is an 
integrable model admitting stable soliton solutions [9]. Experimental observations of 
dark [10, 11] and bright [12, 13, 14, 15] solitons have led to considerable interest in 
understanding the nature of stable collective excitations of BEC in different spatial 
dimensions. A systematic study of the decay of dark soliton has also been carried out 
recently [16]. The above mentioned grey solitons are produced through collision of two 
BECs, when the collisional energy is reduced to a level, where it is comparable with the 
interaction energy and the wavelength of the laser responsible for density modulation is 
larger than the healing length. In this scenario, nonlinearity plays a stronger role than 
dispersion. 

In 1963, in a second quantized formulation, Lieb discovered a collective excitation 
of the condensed bosons, which exhibited a periodic dispersion [17]. The same was later 
identified with a complex soliton at the mean field level [18]. As compared to the dark 
soliton, which necessarily passes through the normal phase, with zero condensate density, 
the Lieb mode can attain the asymptotic density, without reaching the local potential 
maximum corresponding to the normal phase. This is possible because of the complex 
envelope of the Lieb mode, analogous to the Bloch solitons in condensed matter systems. 
Interestingly this non-linear mode has a periodic dispersion very different from other 
collective excitations. Laboratory confirmation of this mode has been hampered by the 
fact that in the experimentally accessible low momenta regime, Lieb mode dispersion 
lies below that of the Bogoliubov excitations [7, 8, 19, 20, 21, 22]. The instability of 
BEC for higher momentum values, where these two modes differ significantly excluded 
this domain from observation. In the recent years considerable attention is paid to the 
study of sound modes not only in Bose gases but also in Fermi systems [23, 24], whereas 
a systematic study with the new experimental developments on the solitonic mode is 
still lagging. 

In this letter, we establish a completely analytical approach to achieve grey 
soliton solution in presence of a modulating harmonic trap and Feshbach induced 
scattering length. Further, a detailed study of their energy and momentum is reported. 
Explicit calculation reveals that, harmonic confinement significantly affects the energy 
momentum profile at low momenta. The fact that, the trap does not allow a uniform 
density and the complex soliton's asymptotic behavior matches with the uniform state, 
leads to this strong interaction between grey solitons and the trap geometry. In contrast, 
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the Bogoliubov mode couples weakly with the trap. This opens up the possibility for the 
coherent control and manipulation of these solitons at low energy, through the temporal 
modulation of the scattering length, as well as the trap [25, 26, 27, 28, 29, 30, 31, 32]. 
We demonstrate this explicitly, by obtaining the self similar grey soliton profile in the 
presence of the harmonic confinement, with time dependent scattering length. The 
solitons are necessarily chirped and can be accelerated, compressed or brought to rest. 
The effect of the trap on the energy and momentum of the solitons is analytically 
computed. 



2. Complex Envelop Soliton 

At low temperature, the BEC, confined in an asymmetric trap, V(x,y) — Vq + V± — 
muj^x 2 + y 2 ) /2 + mu 2 z 2 /2, can be made effectively one dimensional, if cuj_ >> u> . The 
mean field GP equation, with U = 4irh 2 a s /m, 

h 2 - 

ihd t v = (-— v 2 + u\^\ 2 + v)q, (l) 

2m 

then allows separation of transverse and longitudinal dynamics, leading to the 
factorization of the order parameter, such that \1/ = f(z,t)g(x,y,a), where a is the 
local density of particles per unit length [19, 33, 34]. In the weak coupling limit, \g\ 2 
has a Gaussian form and f(z,t) satisfies [7, 8], 

ihd t f(z, t) = -(h 2 /2m)d 2 J + hw ± (l + 2a s \f\ 2 )f(z, t). 

The Lieb mode has been identified with the complex envelope soliton solution: |/(-2)| 2 = 
o"o — o"o cos 2 #sech 2 (,2 cos9/Q [7, 8], where ( is the coherence length and 9 is the Mach 
angle. These modes correspond to localized density dips, which asymptotically match 
with the density cr of the uniform BEC. Explicitly, 9 = sin -1 which makes it evident 



that the soliton velocity u is bounded by the sound speed c s , where c s = \J 2hulA r ^ oas . It 
is also clear that, the shallower solitons move faster. At 9 = 0, where u — 0, the grey 
soliton attains its maximum depth; at 9 = ir/2, u = c s , the maximum attainable speed, 
where a(z,t) = a , corresponding to the uniform BEC. The energy of this extended 
object is analytically computable, with respect to the constant background and is given 
by, §£ (1 - %Y /2 , where E = ^2o- ahu ± cr a ± . 

As will be seen, the presence of the trap strongly influences the grey solitons. This 
arises from the fact that, harmonic confinement does not allow an uniform density and 
the complex soliton's asymptotic behavior matches with the uniform state. This effect 
is particularly prominent in the low momenta regime, where the Bogoliubov modes are 
not significantly affected. In the following, we exactly obtain the self similar grey soliton 
profile in the presence of the harmonic confinement, a time dependent scattering length 
and a phenomenological loss term. This transparently reveals the grey soliton's response 
to the temporal changes of the trap and other parameters. The appropriately scaled 
GP equation, in dimensionless units, can be written in the form [28, 29], 

m = ( - \d zz + 7 (t)H 2 + l -M{t)z 2 + i 9 -^ - (2) 
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Figure 1. Phase step of ir in grey soliton. The parameter values are ctq = 10 7 , 
0.1, A) = 0.5,6» = tt/16. 



K = 



Here, the interaction strength and the spring constant are 7(t) = 2a s (t)/a,B, M(t) = 
Wo(i)/^i, respectively. Further a± = (S/mui) 1 ' 2 and Ob is the Bohr radius. 

The ansatz solution is taken in the form, ip(z,t) = U(z,t) exp[i&(z, £)], where, 



with T = A(t)(z — l(t)) and G(t) = f g(t')dt' . The phase has a quadratic form exhibiting 
chirping: $(z,t) = a(t) + b{t)z - \c{t)z 2 . 

The consistency conditions lead to a Riccati type equation, c t — c 2 {t) = M{t) [29]. 
The width A(t) is given by A(t) = Aoexp(foc(t?)dtf); further j(t) = -f e~ G ' A{t) / A , 
b(t) = A(t) and v(t) = 2A 2 /j,. The center of mass motion is governed by, 



It can be controlled through the trap, condensate motion and scattering length. The 
amplitude is related to the width as, B(t) = ^ A{t). Explicit calculation shows, 
a(t) = do — Jo A 2 (t')dt' . The current conservation yields, = u(l — ^), with 
boundary condition, x' ~ > for a — > o~o, where <7o is the equilibrium density of atoms 
in the moving frame T. Here ft = —2ko§ = \x + A. fi, A are the chemical potential and 
constant parameter controlling the energy of excitation [35] and k = J^. In this frame, 
the density equation can be cast in the convenient form [7], 



U(z,t) = B(t)JoiT)exv[ix(T) + G(t)/2] 



l t + cl - b = Au = V. 




(3) 



The solution for density takes the self similar form [36] 




(4) 




Figure 2. Decrease of soliton depth with increasing velocity in a harmonic trap, 
calculated for 9 = 0,7r/8, 7r/4 and 7r/2, where 6* = sin -1 ^, V being the soliton 
velocity and C s the speed of sound. 



where the Mach angle in the new frame is given by, 

n ■ -1 V ■ -1 U /r\ 

9 = sm — = sm — . (5) 

In the above analysis the temporal variation of the background modifies the sound 
velocity: C s = A^/koq = Ac s [37]. 

The phase step is a direct attribute of the complex soliton, the obtained theoretical 
value is 7r [5], whereas the experimental observation yields 0.77T [1]. The phase change 
across the soliton profile as a function of space-time is depicted in Fig. (1). 

As in the case of the grey soliton without a trap, in the present case also, the 
shallower solitons move faster than the deeper ones, as depicted in Fig. (2). Higher 9 
values correspond to higher speed as is evident from Eq.(5) and correspond to shallower 
grey solitons. 



3. Energy Momentum Dispersion 

The energy of soliton per unit length can be evaluated by subtracting the background 
energy such that E = W — W , with, 

W= J \\(%%) + \l{t)m?) + \M{t)z^]dz 

and ipo = /9 e** < - 2, *' )+G//2 , which corresponds to the BEC with a trivial phase in a trap of 
similar parameter values. The computed energy, 

B = e^V W*-(c* + M)[i^ + 

l 2 (a cos 0] + 2Abua ( cos 6 + b 2 a ( cos 9), (6) 

shows the effect of loss in the exponential prefactor e G . The first term represents 
purely solitonic energy. The term thereafter accounts for the oscillator and chirping 
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Figure 3. Dispersion relation for soliton (blue dashed line) and sound mode (pink 
dashed line), (a) The behavior of Lieb and Bogoliubov mode is shown with 2ir 
periodicity at t = 0. The inset depicts the dispersion relation without the trap at 
t 0. (b) The Lieb mode dispersion for 6 = inpresence of harmonic oscillator trap. 
In the inset depicts the mode behaviour when b ^= 0. 



contributions, which is singular at 9 = n/2. This singularity arises from the physical 
consideration that, the trap does not favour uniform density distribution which excludes 
the asymptotic grey solitons, where a(z, t) = er . The third term, linear in b(t), arises 
from the coupling of BEC momentum with soliton phase and the last term being 
quadratic in b(t), represents BEC translational energy. Fig. 3 (a) depicts the well known 
2-7T periodicity of the Lieb mode in absence of trap. In the inset it has been shown that 
if t ^ 0, the periodicity will be different. The energy and momentum is normalized by 
E and Po, where E = k((A q o~o) 2 and Pq = A c s ao(. The strong effect of the trap at 
low momenta on the grey soliton's dispersion is shown in Fig. 3(b). The presence of the 
translational motion of BEC affects soliton significantly and one observes that in certain 
cases at low momenta, the soliton's momentum can become negative, as depicted in the 
inset of Fig.3(b). There is a point at which the energy is non zero at zero momentum, 
which reveals the rest mass of the solitons. 
The canonical momentum is given by, 




= e G {C s (a (n— - sin 29 - 29) - 2b(a cos 9]. (7) 

The first term can be attributed to soliton, and the second term arises from BEC 
momentum. It is worth mentioning that the soliton velocity can also be computed from 
the hydrodynamic relation: |^ = Au, which matches with the earlier obtained result. 
In the limiting case, i.e., when the trap is switched off, all the expressions match with 
the known results [7, 8]. 

In order to compare the behavior of the Lieb mode and the second sound, we apply 
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a small perturbation [38]: 

8ij)(z,t) = ^/A(t)5p(z,t) exp[i$(z,t)+G(t)/2], 

to a constant background ip®. The modified equation of motion for 5p is given by, 
id t 5p = — \d zz 5p + K,A 2 [2\p \ 2 5p + p^5p*] — i& z 5p. This can be written as an eigenvalue 
equation, Jy t = Hy, where 



and 

H 



-\d zz + KA 2 a -§ z d z 



2' 

Cf> f) _ , 

2 

Assuming a plane wave perturbation i.e., 1/(2, t) = exp[i(/cz— o>£)], the dispersion relation 
can be obtained: 



u = ] j^ + C 2 k 2 + {b- cz)k. (8) 

The Bogoliubov dispersion carries the effect of BEC momentum and chirping in 
the second and third terms respectively. The translational motion of the BEC modifies 
the Bogoliubov dispersion, as has been observed earlier in the case of BEC flowing 
with subsonic [23] and supersonic velocity [23, 39]. In the weak coupling regime, due 
to the presence of the trap, chirping is unavoidable. In the strong coupling scenario, 
Thomas-Fermi approximation accounts for the trap. 



4. Conclusion 



In conclusion, it is observed that, there is a nice and elegant way to extract analytical 
grey soliton solution from a rather complicated and more physically suited representation 
of the Gross-Pitaevskii equation. Explicit calculations show that the harmonic 
confinement significantly affects the energy momentum profile of the grey solitons at 
low momenta, since uniform density is not permitted in the trap. The Bogoliubov mode 
is weakly coupled with the trap. This opens up the possibility for the coherent control 
and manipulation of these solitons at low energy, through the temporal modulation 
of the scattering length, as well as the trap. Explicit analysis of the center of mass 
motion shows that, the soliton can be compressed and accelerated through the trap and 
scattering length. We hope that, the strong coupling of the Lieb mode with the trap 
can be realised through the present experimental setups. 
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